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In this paper we present a generalization of a simple solid-on-solid epitaxial model of thin 
films growth, when surface morphology anisotropy is provoked by anisotropy in model 
control parameters: binding energy and/or diffusion barrier. The anisotropy is discussed 
in terms of the height-height correlation function. It was experimentally confirmed that 
the difference in diffusion barriers yields anisotropy in morphology of the surface. We 
got antisymmetric correlations in the two in-plane directions for antisymmetric binding. 
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1. Introduction 

The surface growth phenomena happen very often, not only in physics. It is the 
case of crystal growthJdQ epidemics, live cell growth, corrosion of materials, the 
spreading of forest fires, eic.cllj Anisotropic growth case could be considered for 
example of the fire as result of influence of strong wind in one direction, or crystals 
growing preferentially by diffusion in an easy diffusion, etc. T n th is paper we would 
like to present some suitable modifications of our earlier modeltB for the anisotropic 
growth case when model control parameters are different in different directions. 

The rules of growth are based on deposition model with particles diffusion limited 
to one-step, to the nearest neighbor onlyB The probability of a particle settling at a 
position is given by the Boltzmann factor exp(—E/kT) where E is particle energy 
at the position, k is the Boltzmann constant, and T denotes temperature. After 
this one inter-atomic-separation long step particle sticks there and stays for rest 
of simulation. We consider three dominant contributions to the particle energy: 
binding to underlying atom energy S < 0, nearest neighbor interaction energy J, 
and diffusion barrier V > 0. We assume that the model control parameters may be 
different in different directions (S x ^ S y or J x ^ J y or V x ^ V y ), which is essential 
to generate the claimed surface morphology anisotropy. 

The nearest neighbor energy J acting between neighbors on same film layer, 
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Figure 1 : Model control parameters S and J. Indexes x and y correspond to x- or 
y-directions of atom pairs. Anisotropy results from J x ^ J y or S x ^ S y . 



had to be enriched by the next nearest neighbor interaction S with underlying 
atoms, see Fig. [j]. (Note that J-term to the atom just below is irrelevant since it 
produces a constant contribution to energy, eventually absorbed by normalization 
factor when probabilities are calculated.) The presence of the S'-term, known as 
the Ehrlich-Schwoebel barrier ,Q'El is necessary to account for the so called edge effect 
when particles do not hop down a terrace edge. 

As we mentioned, the anisotropy will be discussed in terms of the height-height 
correlation function G(s) of a planar vector s, a given distance between the pairs 
of columns 



G(s) = (h(f) ■ h(f+ I)) - (Mr)) 2 , 



(1) 



where (...) denotes spatial average (over all sites r) and (h(r)) is the average height 
of the film. Such a two-site characteristic is necessary for description of anisotropy. 
Note, that G(0,0) = ex 2 , and a is the standard deviation of surface heights. We 
also define dimensionless functions g x = G(l, 0)/G(0, 0) and g y = G(0, 1)/G(0, 0). 
For positive values g x > (or g y > 0), we have terraces- like smooth surfaces in 
x-direction (or y-direction) . On the other hand, for negative values of correlation 
function g x < (or g y < 0), a spiky and rough structures are observed. The 
ellipticity parameter e = g x — g y may be used as a measure of surface anisotropy. 
The ellipticity e vanishes for isotropic surface, anisotropy manifests itself as e =/= 0. 



2. Results of Simulation and Discussion 

Simulations were carried out on 500 x 500 large square lattices with periodic bound- 
ary conditions to minimize boundary effects. The average coverage was ten mono- 
layers (ML). As we mentioned earlier, the anisotropic growth could be provoked 
either by anisotropy in in-plane interaction energy J x ^ J y (see Table and |), m 
binding to underlying layer energy S x ^ S y (Table |3|), and/or by anisotropy in the 
diffusion barrier V x ^ V y when J ^ (Table ||). 

For antisymmetric values of the in- plane interaction energy J x = —J y , the corre- 
lation function in the two perpendicular directions is also antisymmetric, g x = — g y 
(see Table §. In such case particles repel each other in one direction while they 



A simple solid-on-solid model of epitaxial films growth: surface morphology anisotropy 3 

Table 1: Ellipticity e for J x — J, J y — (asymmetric case) for different values of 
diffusion barrier V x — V v — V . S x — S v = 0. 
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Table 2: Ellipticity e for J x = —J y = J (antisymmetric case) for different values of 
diffusion barrier V x — V v — V. S x = S v = 0. 
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attract in the perpendicular one. The surface profile is then similar to the sketch 
presented in Fig. |^c. The anisotropy parameter e does not depend on surface average 
coverage (h) (except of the earliest stages of simulation) as presented in Table fj]. 

The anisotropy parameter e depends on the ratio \J X — J y \/\J X + J y \ of the 
parameters and not on the differences. For the asymmetric case we found that the 
ellipticity e tends to zero for large values of binding energy although J x — J y ^ (see 
Table Hand ^) . The anisotropy vanishes also for relatively large diffusion barrier V 
(Table |l|, | and |). 

Anisotropic growth of the surfaces follows anisotropy in either the atom pairs 
binding energy, or in different diffusion barriers and mobility of atoms in different 
directions. Surface roughness, which is usually characterized by the width of the 
surface, is not a suitable choice for anisotropy discussion as a single site property 
only. Instead we need a two-site properties such as the height-height correlation 
function G(n,m) for columns separated by a vector (n,m), or ellipticity parameter 
e. The main results of the simulations may be summarized as follows: 

Firstly, for the isotropic case the correlation function G(1,0) or G(0, 1) itself 
reflects the roughness. For random deposition G(1,0) = G(0, 1) = 0. For flatter 
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Table 3: Ellipticity e for S x — S, S y = for different values of diffusion barrier 
V x = V y = V. J x = J y = 0. 
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Table 4: Ellipticity e for V x = V, V y = for different values of binding energy 

Jx Jy J • S x Sy 0. 
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Roughness a [ML] and ellipticity e for asymmetric anisotropic case (J x = J 
--J-10). 

~1 TToo Too Too Too ooo - 

a 0.95 0.94 0.99 1.07 1.78 
e -0.08 -0.09 -0.16 -0.28 -1.03 



Roughness a [ML] and ellipticity e for asymmetric anisotropic case (J x = J 
--J+10). ' 

J 0.00 1.00 2.00 5.00 10.00 

a 19.04 20.63 21.85 22.55 22.57 

e 0.75 0.31 0.11 0.00 0.00 



Table 7: Dependence of dimensionless correlation functions in perpendicular direc- 
tion (g x ,g y ) for antisymmetric case J x = —J y = —2.0 for different film average 

heights (h) [ML], 

(h) 1 5 10 25 50 75 100 500 
g x 0.54 0.75 0.79 0.80 0.80 0.81 0.80 0.80 
g y -0.45 -0.73 -0.80 -0.81 -0.82 -0.83 -0.83 -0.84 
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Table 8: Correlations expressed in G(0,0) units for anisotropic case (J x = —J y = 

-2.0). 

n 1 2 3 4 5 

G(0,n) -0.80 0.49 -0.24 0.10 -0.03 
G(n,0) 0.79 0.53 0.31 0.15 0.06 
G{n,n) -0.66 0.29 -0.09 0.01 0.00 



Figure 2: Schematic sketch of some typical situation: (a) smooth surface, (b) rough 
surface, (c) anisotropic surface. 



surfaces we claim positive G(1,0) and G(0, 1). This is so since any two adjacent 
sites are expected to have the same heights and so the deviation from the average 
film thickness are of the same sign resulting in the positive contribution to the 
correlation. And on the contrary, rough surfaces show negative G(l, 0) and G(0, 1). 
This, however, does not yield information on possible anisotropy of the surface when 
G(1,0) ^ G(0, 1). Then anisotropy parameter e properly normalized to G(0, 0) is 
used. 

In the limiting case of large diffusion barriers the migration stops and again 
we restore the random deposition model with its consequences, correlations G(s ^ 
0) = and anisotropy e = 0. However, if migration takes place then G(s ^ 0) ^ 
and we may expect that only for distant sites G(n,m) 0, when n — > oo and/or 
in — ► oo. In fact, we often observed a damped character of G(n, 0) dependence 
on the distance n between atoms in the x-direction (or similarly G(0,n) in the y- 
direction). The roughly exponential character of the decaying G(n,0) and G(0,n) 
on n may be observed for the direction which shows a smooth variations in column 
width (negative J, small a). The decrease of G(n,0) and G(0,n) on n may have 
an oscillatory character changing the sign of G(n,0) or G(0,n) if rough and spiky 
structure (positive J, large a) is observed. Therefore anisotropic case J x = —J y 
produces different G(n, 0) and G(0, n) on n dependence in x- and y-directions, 
exponential decrease along x-axis for atom pairs correlation G(n,0), and damped 
oscillatory dependence of G(0, n) in y-direction, as it is shown in Table ||. 

The close-distance correlations G(l, 0), G(0, 1) and G(l, 1) may also be helpful 
to distinguish between different classes of surface morphology. As we mentioned 
earlier, G(1,0) < (G(0, 1) < 0) indicates rough surface, G(1,0) > (G(0,1) > 0) 
indicates flat surface. Therefore, we expect for an isotropic flat surface as in Fig. ||a 
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(7(1,0) = G(0, 1) > and G(l, 1) > 0. For a chessboard-like surface, Fig. |b, we 
get G(1,0) = G(0, 1) < and again G(l, 1) > 0. Thus G(l, 1) cannot distinguish 
between isotropic smooth (Fig. ^a) and isotropic rough (Fig. |^b) case. However, for 
anisotropic case, Fig. ||c, of a terrace-like structure we get G(l, 0) > 0, G(0, 1) < 
and G(l, 1) < 0. We conclude that combination of all G(l, 0), G(0, 1) and G(l, 1) 
information is indicative of the type of the film morphologies, its roughness and 
anisotropy. 

The anisotropy in surface morphology, particles migrations and accommodation 
coefficients, as well as surface magnetic properties, were verified experimentally. 
For instance, Mo et al worksi^Ej show a highly anisotropic diffusion of Si and 
Ge atoms on a Si(001) surface, with the easy diffusion direction along the dimcr 
rows. Surface migration of atoms is at least one thousand times faster along the 
substrate dimer rows than in the direction perpendicular to thcm.BIll It could be 
explained theoretically by assuming^ the difference in energy activation for atom 
diffusion AE = E\\ - E x « 0.4 eV.El 
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